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ABSTRACT 

The problem of the pairing interaction is dealt with even Grass- 
mann variables in the hamiltonian framework . Eigenf unctions of 
given energy, seniority and zero third component of angular mo- 
mentum are given in terms of single particle and collective bosonic 
fields. 

Si tratta il problema dell'interazione di pairing con le varia- 
bili di Grassman pari nel formalismo hamiltoniano . Si danno le 
autofunzioni di data energia, senioritd e terza componente del 
momento angolare pari a zero, in termini di campi bosonici di 
particella singola e collettivi. 



Recently the problem of bosonizing the pairing hamiltonian Hp has been 
addressed in the formalism of the path integrals over Grassmann variables 
[|I|]. In this context, using even elements of the Grassman algebra, a func- 
tional integral representation of the correlation functions has been given by 
performing a change of variables in the Berezin integral. It has thus been 
possible to express the pairing action and the ground state wave function 
through "collective variables", namely specific linear combinations of the 
even elements of the algebra, and to obtain the ground state energy of Hp. 

However, to solve the problem of the excited states of Hp (with non zero 
seniority v) in the path integral framework has proved to be a quite difficult 
task to perform. Accordingly here we address the v 7^ problem resorting 
to the Grassmann algebra in the hamiltonian framework. 

Of course the spectrum of Hp has already been found long time ago within 
the quasi spin formalism , but still we believe it useful to address the same 
problem in this new context to shed light on the role of composite fields in 
nuclear physics and to provide explicit expressions for the eigenfunctions of 
Hp with a non vanishing seniority. 

To start with, we shortly recall the essential elements of the hamilto- 
nian many-fermions problem in terms of Grassmann variables (the bosonic 
problem, likewise, can be formulated in terms of holomorphic variables) || . 

One exploits the isomorphism between the Fock space F generated by the 
fermionic creation operators af, • • • a^" n and the algebra Q + generated by the 
set of totally anticommuting objects A*, ■ ■■X^n (odd elements of the algebra). 
The isomorphism is defined by mapping the vectors df...dj~\0 > onto the 
elements A*...A*j. The image of a generic vector >G F under this mapping 
will be denoted by \I/(A*). 

The scalar product in F in terms of Grassmann variables reads then 

< *!|* a >= J dA* n dA 2a ...dA*dA 1 (^ 1 (A*))>_(A*A)* 2 (A*) ) (1) 
the measure being 

/i ± (A*A) = e ± ^ A * A \ (2) 

Actually to define the kernels of integral operators the measure with the plus 
sign is needed. 

Indeed a linear operator in normal form, namely 
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0= E E :, - h -' - JI "n-","r-" l , , (3) 

can be expressed in terms of Grassmann variables according to 

0(X*,X)= E ^ 1 --"- r - A A:...A: i A,,..^ , (4) 

h,...ik,ji---3k 

and its associated kernel reads 

K (A*,A) = 0(A*,A) /U+ (A*A). (5) 
Then the action of O on a state \1> is given, in integral form, as 

(0*)(A*) = y [rfA , *rfA , ]K (A*,A , ) / u_(A , *A')^(A'*) . (6) 

We apply the formalism, shortly revisited above, to a system made of one 
kind of nucleons (e.g. neutrons) in a level specified by the angular momentum 
j, hence the indices labelling the A variables will correspond to the third 
components m of j. 

Now if only pairs of nucleons coupled to a total angular momentum with 
J z = are considered, then the exponentials appearing in fi± can be trans- 
formed as follows 

j [d\*\ exp (±A^A m ± y m \- m ) 7(0*) = / W] exp (C0 m ) /(</>*) , (7) 
in terms of the even Grassmann variables 

<fim — A_ m A m . (8) 

For N (even) nucleons in a single particle level j, the pairing hamiltonian 

reads 

H$ = -gpA^Aj , (9) 

where 

[tt j j 
A i = \~o ^ (jm,j -m\00)a j _ m a jm = E m %-m%m , (10) 

V m=-j m=l/2 
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being Q = (2j + l)/2 and (jm, j— m|00) the usual Clebsch Gordan coefficient. 
In the frame of the Grassmann algebra, we recast instead Eq. (§) as follows 



Hp = -g P A*A 3 (11) 

with 

*i = E (-l) J " m A- m A m = (-l) J_m 0- • (12) 

m=l/2 m=l/2 

We explore now whether the 0's can be the building blocks of the eigen- 
states of Hp\ Indeed this turns out to be the case, since, according to @ 
and (0), the action of the hamiltonian over states set up with the 0's is 

= [Wdr])^ (</>*, r/)e£ ,f 'V(7f) , (13) 

the kernel reading 

Kf(4>*, r?) = i#(<f , r?)eS^ . (14) 
Hence the eigenvalues equation 

H^i>{(f>*) = J [drfdr^Kf (</>*, rfieE^tf) = £ (j 'V(0*) (15) 

follows. 

We then look for eigenstates of ifp^ in the form of linear combinations 
of products of 0's. More specifically, for fixed Q and for n = N/2 pairs, the 
wave function will be a superposition of (^j monomials, of the type 

(2) 

#n=E»m 1 "--CJ i , (16) 

8=1 

where the are complex coefficients. 

In order to study the secular equation of Hp\ we construct the matrix 
associated to the latter in the basis of the states flTE|). To start with, we 
consider one pair only. In this case, the wave function and the eigenvalue 
equation read, respectively, 

i>m = E (-ir m A.c (17) 

m=l/2 
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and 

(8 + l)(3 m + £ Pp = ° ( 18 ) 

p(^m)=l/2 

where £ = E/gp . Moreover the representation of Hp in the basis (|T7| ) is 



/ 8 + 1 1 1 

1 8 + 1 1 

1 1 8 + 1 

1 1 ■■■8 + 1 1 

1 1 8 + 1 



V 1 



1 



1 



1 \ 
1 
1 
1 
1 

8 + 1 J 



(19) 



a matrix of dimension Q with the index m labelling the rows and the columns. 
This matrix is obviously invariant for a permutation of the values of m. 
Elementary methods lead then to the characteristic equation 



{8 + fi) 8 



n-i 



. 



We thus see that, of the Q expected real roots, two only are distinct, namely 
the lower one 

S ( n =l) = -n 



with multiplicity 5=1 and 



8 2 {n = 1) = 



with multiplicity 5 = Q — 1. 

In turn, the associated eigenvectors read, respectively, 



j-m,* 
m 



and 



j'-i 



m=l/2 



1 



(20) 



(21) 



* 2 (n=l)=JV £ /? m 

m=l/2 

A/" being a normalization factor. The above findings agree with the physics 
of the hamiltonian Hp\ which acts only in presence of J = pairs. Also the 
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eigenvectors (^(], ^1]) are orthogonal and the eigenvalues are those obtained 
in the usual quasi-spin framework. 

Of significance is that the characteristic equation for the case with n = 
Q — 1, namely 

(£ + Q-l)(3 m + £ A> = ' ( 22 ) 

p(^ m )=l/2 

has the same structure of (13) ■ However the states read now 



W) = E (-l) i_m An«/2 

771=1/2 



^m-lVm+l 



the index m corresponding to the missing pair. 

By extension, one finds that the case with n pairs present is equivalent 
to the one with Q — n pairs: indeed it leads to the same matrix but for 
the diagonal elements. Hence in the following we shall confine ourselves to 
explore the cases with n < ^ only. 

Next we address the case n = 2. Here the wave function reads 

j m— 1 



entailing the secular equation 



E E( 

771=1/2 71=1/2 



\3-rni 



(23) 



(£ + 2)(3 mn + Yl (fa™ + Ppn) = • 
p(^ra,m)=l/2 



(24) 



As for the case n — 1, the associated matrix can be written in ( ! equivalent 
ways (not all of them being distinct), each one corresponding to a different 
labelling of the states, but yielding the same determinant. 

For example, for Q = 4 (the smallest Q for n = 2) one finds the following 
6x6 matrix 



(8+2 1 1 1 
1 £+2 1 1 

1 1 8 + 2 
1 1 £+2 
10 11 
V 1 1 1 



1 \ 

1 

1 1 
1 1 

£ + 2 1 

1 £ + 2 j 



(25) 
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with the characteristic equation 

(£ + 6)(£ + 2) 3 £ 2 = . 

Thus, in this case, out of 6 real roots only 3 turn out to be distinct, namely 
the lowest 

£ (n = 2) = -6 
with degeneracy 5=1, the intermediate 

£ 2 {n = 2) = -2 

with degeneracy 5 = 3 and the highest 

£ 4 (n = 2) = 

with degeneracy 5 = 2. The associated orthogonal eigenvectors are found to 
be 



(26) 



m Q (n = 2) = 27V (-01/203/2 + 01/205/2 - 01/207/2 
-03/205/2 + 03/207/2 ~ 05/207^) 

m 2 (n = 2)=M 2 [ax (0t/ 2 03/ 2 - 05/207/ 2 ) 
+a 2 (0 1/2 05/2 - 03/207/2) 
+a 3 (0^/207/2 - 03/205/2) _ 

m A (n = 2)= Ma [(62 " 63) (01/203/2 + 05/207/ 2 ) 

+b 2 



(27) 



* J* 1 J. * J* 



+ 



+h (01/207/2 + 03/205/2 



(28) 



Ao, A2, A/4 being normalization factors and a\, a 2 , a^, b 2 and 6 3 free parame- 
ters. 

Of significance is the correspondence between the number of parameters 
in the wave function (in eq. (|26|) the parameter has been absorbed in the 
normalization factor) and the roots degeneracy, namely the number of inde- 
pendent eigenf unctions of given energy, seniority and zero third component of 
the angular momentum which can be set up with the 0's as building blocks. 
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A few comments are now in order. First, clearly, three distinct eigenvalues 
are just expected for two pairs, corresponding to zero, one and two pairs 
broken respectively. By extension, no matter what the value of Q is, the 
distinct roots of the secular equation will always be n+1. Their degeneracy 
can be found, as above stated, by inspecting the associated wave functions. 
However it can be understood on the basis of the following counting rule: 
the lowest eigenvalue is always unique 



£ (n) 



6n 



1 



and corresponds to the most collective eigenvector \l/ ( see Eqsj20|-26). The 
second eigenvector ^2 has a broken pair: since the states available for a pair 
are Q, the broken pair can sit in any of them, except for the one corresponding 
to ^o- hence its degeneracy is 



S 2 (n) 



So 



1 = n- 1 



The next eigenvector \l/4 has two pairs broken. The associated degeneracy is 
obtained by subtracting, out of the states available for two pairs, the states 
already occupied by \l/o and ^2, namely 



SJn) 



64 



1 - 



1 



In general the following degeneracy for the eigenvalues corresponding to 
seniority v is found to hold (of course, a binomial coefficient with a negative 
lower index vanishes) 



' n N 
J- 2, 



SJn) 



n 



(29) 



v being an even non-negative number referred to as seniority, to comply with 
the existing literature. 
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Finally we explore the general case. Here the structure of the symmetric 
matrix of dimension associated with Hp^ is easily found to be 



f £ + n V 1 \ 



^ V 1 £ + n J 



(30) 



where the symbol V 1 means that the upper (and the lower) triangle of the 
matrix is filled with zeros and ones. Indeed the matrix elements of Hp^ turn 
out to be one, when the bra and the ket differ by the quantum state of one 
(out of n) pair, otherwise they vanish. The diagonal matrix elements simply 
count the number of pairs. 

An elementary combinatorial analysis shows that the number of ones in 
each row (column) of the matrix is given by n(Q — n). Indeed a non vanishing 
matrix element has the row specified by n indices whereas, of the indices 
identifying the column, n — 1 should be extracted from the n ones fixing the 
row in all the possible ways, which amounts to n possibilities. The missing 
index should then be selected among the remaining Q — n ones: hence the 
formula n(Q — n) for the number of ones. 

No matter where these are located Q, in the determinant of the matrix 
the first row (or column) can be replaced with the sum of all the rows (or 
columns). This leads to a new determinant with, e.g., the first row (or 
column) filled with ones and having the expression £+n+n(Q— n) factorized. 
Hence 

1. the lowest eigenvalue is given by the integer 

£ = - n (i + n- n ) = -n-n(fi-n) ; (31) 



2. it is fixed uniquely by n and by the number of ones present in each row 
(or column) of the matrix. 

Next notice that, setting £ = in (|30"D , for a given Q, the number of pairs 
n(< ^), beyond fixing the dimension of the matrix and the number of ones 

* As already noticed, any ordering of the states should of course lead to the same results. 
There are ( ) ! ordering alltogheter. Note however that for ft and n large the number of 
matrix with n(Cl — n) ones in each row and column is larger than (^)I- 
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in each row (column), also represents the common value of the diagonal 
elements. We have numerically checked that the determinant of such a matrix 
always vanishes. Hence it follows that 8 = is always an eigenvalue of Hp^ 
and actually the upper one. Indeed it corresponds to the situation where all 
the pairs are broken. 

Let us now insert into the matrix flSCl) the lowest eigenvalue (PI). We 
obtain 

/ —n(Q — n) 



det 



) the lowest eigenvalue 
V 1 \ 



V 







(32) 



V 1 -n(tt -n) ) 

We extensively checked by direct computation that the above determinant 
continues to vanish under the replacement 



n) 



—n(Q — n) + 



n 



+ 1 



2 V 2 

Hence the characteristic equation in the general case turns out to be 



n 

v/2=0 



£ + n(Q-n+l] 







(33) 



and the whole spectrum of the pairing hamiltonian follows. 

The associated wave functions will be linear combinations of r^j mono- 
mials, each monomial being the product of n 0's, as usual. 

In conclusion, we like to shortly address the problem of the angular mo- 
mentum of our eigenstates, which belong to a definite value of energy, senior- 
ity and third component of the angular momentum, but not of the angular 
momentum. However, their building blocks, the monomials, can be expressed 
as superpositions of "collective" bosons $j of definite angular momentum, 
according to 

( j> m = V2j2(jm,3-m\J0)$ J . (34) 
j 

It should be realized however that the nihilpotency of the variables 0's in- 
duces a set of constraints to be fulfilled by the "collective" variables $j. As 
an example, consider again the case Q = 4, n = 2. Here the four constraints 
are translated in the following four equations for the variables $j 



i> 2 
in 



with J 



0,2,4 and 6: 



$1 



(35) 
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77V2T$o$2 = -88$^ - 10$^ + 98$§ - 44V33$ 2 $ 4 - 35V21$4$6 (36) 
231v / 77$ $4 = -726$ 2 + 432$ 2 + 294$ 2 - 

20V33$ 2 $4 - 105v / 77$ 2 $6 + 280v / 21$ 4 $6 (37) 
llv^^e = 20$^ - 20$ 2 1 - 5v / 33$ 2 $4 - 4 v / 77$ 2 $ 6 + 4 V / 21$ 4 $6 . (38) 

Hence the v = eigenstate ([ 



when expressed through the $j, reads 



# (n 



No- (-3$ 



(39) 



which, exploiting (p5|), actually becomes 



^o(n = 2) = -2A/" $2 . 



(40) 



The above is the product of two s-bosons, as the variable $o is often referred 
to, or of the wavefunctions of two unbroken pairs. 

Analogously the v = 2 states, exploiting the constraints (|3"6"|-|5E|) , turns 
out to read 



N 2 



V3\/7\/ll 
+ (-3ai - 2a 2 - 8a 3 ) V3%<5> 



(4oi - 2a 2 - a 3 ) Vll$o$2+ 

(-2 ai - 5a 2 + 2a 3 ) V^o^el , (41) 



which shows that a particular choice of the parameters allows to select a state 
of good angular momentum (J = 2, 4, 6). Hence this state can be reduced to 
the product, e.g., of an s and d boson, the latter corresponding to the wave 
function of the broken pair. 

More subtle is the v = 4 case. Here the procedure above outlined leads 
to the following expression of the wavefunction 



Na 



h f --$2 - — $4 + — $fi + 

V 7 2 77 4 11 6 7V33 



+63 



7 



7 



40 



7V33 



$ 2 $ 4 



$2$4 + ^p^2$6 + 

4 4 \ 

$ 2 $ 6 + ^=$ 4 $ 6 
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llV2T 



$4$f 



77 



(42) 



where, remarkably, the collective $0 does not appear, as it should, since 
is associated with two broken pairs. 
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